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ABSTRACT 

We determine the characteristic length scale, L(t), in phase ordering kinetics for 
both scalar and vector fields, with either short- or long-range interactions, and with 
or without conservation laws. We obtain L(t) consistently by comparing the global 
rate of energy change to the energy dissipation from the local evolution of the order 
parameter. We derive growth laws for 0(n) models, and our results can be applied 
to other systems with similar defect structures. 
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Systems quenched from a disordered phase into an ordered phase do not order 
instantaneously. Instead, the length scale of ordered regions grows at a character- 
istic rate as different broken symmetry phases compete to select the ground state. 
During this process, energy is dissipated and topological defects, if present, are de- 
stroyed. Traditionally, systems with scalar order parameters, such as binary alloys 
and Ising models, have been studied 01- I n such systems domains of both phases 
grow and the intervening domain walls, which are the characteristic topological de- 
fects of scalar systems, shrink in total area; the energy of the system is dissipated as 
the total area of domain wall decreases. Recently, there has been a growing interest 
in systems with vector and more complex order parameters which have defect struc- 
tures such as lines, points, and textures. An exciting series of experiments 0, |J and 
simulations f| |, [7|, || |] has been performed to explore these systems. Of particular 
interest is the time dependence, L(t), of the characteristic length scale. This growth 
law provides a basis for further analysis and a fundamental test of approximation 



schemes. Growth laws have been obtained in individual cases |fL0|, |TT| , |T|, |13 |, but 
the variety of possible systems has prevented the development of a unified approach. 
Renormalization group (RG) ideas |L2|, |14| have been used to predict the leading 
power-law dependence of growth laws, but are not easy to apply in systems without 
conservation laws. An approach to calculating growth laws, including logarithmic 
factors, for a wide variety of systems would be welcome. 

In this Letter we determine the characteristic length scale, L(t), for quenched 
systems with both scalar and vector fields, with either short- or long-range interac- 
tions, and with or without conservation laws. We base our approach on the energy 
dissipation that occurs as the system relaxes towards its ground state. We evalu- 
ate the energy dissipation by considering topological defects, which dominate the 
dynamics when they exist. We compare the global rate of energy change to the 
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energy dissipation from the local evolution of the order parameter. For systems 
with a single characteristic scale describing correlations as well as dissipation we 
self-consistently determine the growth law, L(t). 

Our results are summarized in tables 1 and 2. We find that conservation laws 
(/i) and long-range interactions (a) are irrelevant to the growth law below and 
above, respectively, values determined by the number of components (n) in the 
order parameter of the system. At the marginal values, logarithmic factors are 
introduced. The growth laws obtained are independent of the spatial dimension 
of the system. Our results should be compared to growth laws obtained by other 
means. Heuristic arguments yield growth laws of L ~ i 1 / 2 for non-conserved JXTj] and 
L ~ t 1//3 for conserved [jnj scalar order parameters. RG arguments [T2[] reproduce 
the t 1 / 3 law for conserved scalar fields, predict L ~ t 1 / 4 for a conserved vector field, 
and also treat conserved systems with long-range interactions |L4]. Simulations f| 
suggest a L ~ t 1 ^ 2 growth law for non-conserved vector fields, and simulation || and 
experiments || obtain the same law for the related tensor fields that describe nematic 
liquid crystals. This result has also been obtained by approximate treatments for 
general n jlS] and verified in the large- n limit [13|]. Our approach provides a unified 
theoretical basis for the leading power-law behavior suggested by these results, and 
also determines any logarithmic factors. 

The dynamic scaling hypothesis describes most systems in the late stages of 
growth |TJ|. Accordingly, the correlation function of the order parameter, C(r,t) = 
^0(x, t) ■ 0(x + r, t)^, with an average over initial conditions, should exhibit the scal- 
ing form C(r, t) = f(r/L(t)), with a single characteristic length scale L(t). Fourier 
transforming this, we obtain the scaling form for the structure factor 5(k, t) = 

4w • 0-k(*)), 

S(k,t)=L(t) d g(kL(t)) , (1) 
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where d is the spatial dimension. A generic energy functional of such a system, with 
an n-component field, <p(x), is 

H[$\ = J d d x [(V0) 2 + V{4>)\ , (2) 

where V(4>) is a "mexican hat" shaped potential such as V(4>) = (4> 2 — l) 2 - The 
equation of motion for the ordering kinetics of the Fourier components 0k is 

dt^-kHdH/d}-*), (3) 

where we work at temperature T = 0, with no thermal noise, since T appears to be 
an 'irrelevant variable' for ordering kinetics within the ordered phase fl2], [13]]. The 
conventional non-conserved (model A dynamics) and conserved (model B dynamics) 
cases are \i = and \i = 2, respectively. 

Integrating the rate of energy dissipation from each Fourier mode, and then using 
the equation of motion we find 

de/dt = ^({dH/dfo-dtfa) 

= ~J k k ^ (dt$k-dt?-*) , (4) 

where e = (H) /V is the mean energy density, and / k is the momentum integral 
/ d d k/(27i) d . We will relate the scaling behavior of both sides of @ to that of 
appropriate integrals over the structure factor, S(k, t). Either the integrals converge 
and the dependence on the scale L(t) can be extracted using the scaling form (P, or 
a divergent contribution dominates the integral. Any divergence will occur because 
of small-scale structure at kL ^> 1, since no structure has formed yet for kL 1. 
This small-scale structure will only come from topological defects, so we expect 
S(k,t) to be proportional to the density of defect core, pdef ~ L d ~ n / L d ~ L~ n for 
kL ^> 1 [IB . The scaling form (|l|) then implies 



S(k, t) ~ L- n k~ {d+n) , kL > 1, (5) 



which is a generalized Porod's law for n < d [HI, ED . 



We first calculate the scaling behavior of the energy density, e, which is captured 
by that of the gradient term in (||]): 



= / k 2 L d g(kL) , (6) 

where we have used the scaling form (P for the structure factor. If the integral 
converges, then a simple change of variables yields e ~ L~ 2 . If the integral diverges 
for kL ^> 1, then we use Porod's law (||) and impose a cutoff at k ~ l/£, where £ is 
the defect core size. We obtain 



e ~ L- n C~ 2 , n < 2 , 

~ L- 2 ln(L/0 , n = 2 , 

~ L~ 2 , n > 2 . (7) 

We see that the energy is dominated by the defect core density pdef for n < 2, by the 
defect field at all length scales for n = 2, and by variations of the order parameter 
at scale L(t) for n > 2. 

We now evaluate the right side of (H) in a similar way. Using the scaling assump- 
tion for the two-time function, ((fk(t) ■ 0_k(t')) = k~ d g(kL(t), kL(t') ), we find 

d 2 



dtdt' 



t=f 

\2 T d~2 



= (dL/dtyL d - 2 h(kL). (8) 

When the momentum integral on the right of (Q) converges, we obtain, using (^), 
de/dt ~ —{dL/dt) 2 L^~ 2 . However, if the integral diverges we need to know the 
behaviour of @ for kL ^> 1. 

In this regime, where the important structure is short compared to the typical 
interdefect spacing, the scaling functions are determined by the field of a single 



defect |TJ|, and the defect field will be comoving with the defect core, 

d t <p = uj x <p — v ■ . (9) 

Here, v is the velocity of the defect and uj represents rotations for vector fields. For 
kL 3> 1, the velocity term dominates because the rotation scales as uj ~ (dL/dt)/L 
while the gradient term involves the short-scale structure, v ■ V ~ (dL/dt)k. In 
estimating v ~ (dL/dt), we are assuming that dissipation occurs by the evolution 
of defects of scale L(t). Using d t (f>k ~ (dL/dt) kifik and Porod's law (|), we have 

(dtfk ■ d t 0- k ) ~ L- n k-( d+n ~V (dL/dt) 2 , kL > 1 . (10) 

Using ( P~D|) in the right of (Q) shows that the integral is convergent for kL ^> 1 when 
n + ji > 2. Otherwise the integral is dominated by k near the upper cut-off l/£. We 
obtain 

J k-» (dtfa ~ L~ n C +fM ~ 2 (dL/dt) 2 , n + /i<2, 

~ L~ n ln(L/0 (dL/dt) 2 , n + /i = 2 , 
~ L^' 2 (dL/dt) 2 , n + /i>2. (11) 

In what cases does dissipation occur primarily at scales as we have just 

assumed? For n > 2, the energy density (|7p, and hence dissipation, is dominated by 
variations at scale L(t). For n < 2, with d > n, the energy density is proportional to 
the defect core volume, but the dissipation is still dominated by scales of order L(t). 
For this case the rate of energy dissipation at a length scale I is given by the rate at 
which the core volume disappears at that scale. So the total energy dissipation is 
given by 

de/dt ~ J dln(l)U d - n -\ (12) 

where I is the rate of evolution of defect structures of scale I, n(l) is the number 
density of features of scale I, and I l d ~ n ~ l is the rate of change of the defect core 



volume at scale I. But ~ n(l) I is just the number flux of defect features, and 
imposing the continuity equation we see that N ~ j(0), where N is the average 
density of defects. N scales as a volume, N ~ ^/L d , and so iV is a constant for 
times of order L/L. This implies that is constant for I <C L, in order to provide 
a constant rate of defect extinction. We see that de/dt ~ / dll d ~ n ~ 1 for I <C L, and 
the convergence at small I for d > n implies that structures with scales I ~ L(t) 
dominate the energy dissipation. 

For n = d with n < 2, the marginal case in the above argument, the energy den- 
sity (0) is dominated by the core energy of point defects and dissipation is dominated 
by defect pairs annihilating. Since the dissipation occurs at separations ! ~ ^ < L 
we do not expect our approach to cover these cases. In fact, since the bulk energy 
density does not depend on defect separations I ^> £, we expect the system to be 
disordered, with an equilibrium density of defects, at any non-zero temperature. At 
T = we expect slow growth laws that depend on the details of the potential V (</>). 
These cases, n = d with n < 2, are at their lower critical dimension. 

The 2d XY model is outside the scope of this Letter, since a more detailed 
argument |^TJ shows that dissipation occurs significantly at all length scales between 
£ and L(t). Naively, we might expect scaling violations for this case. Yurke et al. 
suggest a growth law of L ~ (t/ln£) 1//2 for non-conserved order-parameters, while 
simulations by Mondello and Goldenfeld [|7| find L ~ t 1 / 4 for the conserved case. 

In the cases where dissipation occurs at the scale of L(t), we compare the rate of 
energy dissipation between ( [TTj ) and the time derivative of (0) to obtain dL/dt and 
hence L(t). The results are summarized in table 1. 

For non-conserved fields (fj, = 0), we find L ~ t 1 / 2 for all systems (with d > n 
or n > 2). Leading corrections in the n = 2 case are interesting: the InL factors in 
(13) and (17) will in general have different effective cutoffs, of order the core size £. 
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This leads to a logarithmic correction to scaling, L ~ t l l 2 {l + 0(1/ hit)), and may 
account for the smaller exponent (~ 0.45) seen in simulations of 0(2) systems |g, ||. 

For conserved fields (fi > 0) our results agree with the the RG analysis [fL2| , with 
additional logarithmic factors for the marginal cases n = 2 and n + fi = 2. We see 
that the conservation law is only relevant for n + fi > 2. Therefore for vector fields 
(n > 2), any > is sufficient to change the growth law, while for scalar fields 
(n = 1) the conservation law is irrelevant for /i < 1, in agreement with the RG 
analysis [|12|]. Simulations by Siegert and Rao || for n = 2, with d = 3 and = 2, 
obtain growth exponents slightly over 1/4, which may be difficult to distinguish 
from our predicted L ~ (tint) 1 / 4 behavior. 

The above discussion covers systems with purely short-ranged interactions. We 



can also add long-ranged interactions to our energy functional (H) [H, [14 

,(0(x + r)-0(x)) 2 



H LR ~ I d a x I d a r 



(13) 



The evaluation of the right side of (f|) is unchanged, but the mean energy is now 
given by 

[k ff S(k,t) , a<2. (14) 

We consider o < 2 only, since for a > 2 the energy is dominated by short distances 
and we recover the short-range k 2 form [23]. Inserting the scaling form ([!]) for 
^(k, if:), and noting the asymptotic behavior @, gives 



~ L- ff ln(L/0 
~ L a 



n < a, 



n = a. 



n > a. 



(15) 



We then compare the rate of energy dissipation between ([TT|), which is unchanged, 
and the time derivative of (|15D, to find dL/dt and generate table 2. Our results agree 



with the leading power laws from an RG treatment by Bray for the conserved 
case (/i > 0), while here we obtain the growth laws for all cases, with any logarithmic 
corrections. 

From fll5D , we see that the long-range interactions are relevant only for n > a 
fl4| . This correctly indicates that defects are non-interacting for n < a, where the 
dynamics are driven by a generalized surface tension. Note that cases with a < 
lead to an instability at fixed non-zero wave-vector which breaks the scaling (p]) and 
hence our approach does not directly apply. It is also interesting to note that naive 
power counting using Hlr ([13]) in the dynamics (|3]) yields the incorrect growth law 
L ~ i 1 / " for non-conserved scalar fields. The fact that power counting obtains the 
correct growth law, L ~ t 1//2 , for short range interactions with n < 2, can be seen 
as a fortuitous balancing of the UV divergences in (|) and (H). 

Since systems without topological defects (n > d + 1) will have convergent mo- 
mentum integrals for kL 1, we obtain L ~ t 1 /( fT +^) ; where a = 2 for short-ranged 
forces. We can also apply this result to systems with topological textures {n = d+1), 
even though the appropriate Porod's law is not known. Since defects with n > d 
must be spatially extended and without a core, they will have a smaller large-/c tail 
to their structure factor S(h, t) than any defects with cores. So for n > 2, when 
the energy dissipation clearly occurs at length scales of order L{t) (see (0)) and 
the momentum integrals for defects with cores converge, our results should apply. 
As a result, tables 1 and 2 will apply for any system, except perhaps those with 
d < n < 2 [|24|| . Of course, whether the structure factor scales, as in ([]]), for systems 
with textures is an open question. 

The strength of our approach is that it can be applied to systems with more 
complicated order parameters than n-component vectors. The details of the energy 
functional @ are unimportant to this approach |2|| . All we need is the existence of 
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some short or long-ranged "elastic energy" (a), an effective conservation law (//), and 
the defect structure if any. Our derivation is independent of the initial conditions, 
and so, e.g., applies equally to critical and off-critical quenches as long as the system 
scales at late times. We simply choose a Porod's law (§) to represent the dominant 
defect type, which is the one represented by the smallest n. For example, in bulk 
nematic liquid crystals, the existence of string defects leads to (|5|) with n = 2, which 
with no conservation law implies a L ~ t 1 ! 2 growth law, consistent with recent 
experiments [[| and simulations ||. 

In summary, we obtain growth laws for phase ordering based only on the dynamic 
scaling hypothesis and a generalized Porod's law for the large-momentum tail of the 
structure factor. By focusing on the total energy dissipation, rather than the details 
of the system, we do not need to consider defect dynamics. This results in a powerful 
yet simple method that addresses many cases of interest. 

We thank R. E. Blundell, D. A. Huse, H. H. Lee, and M. A. Moore for discussions. 
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L(t) 


n<2 n = 2 n>2 


n + 11 < 2 
n + /i = 2 
n + /i > 2 


t l/2 

(t/hit) 1 ' 2 t i/2 

t l/(n+ M ) (tlnt) 1 /^) iV(2+M) 



Table 1: The growth law of the length scale L(t) for various number of components, 
n, and conservation laws, /i. Note that n = d is excluded for n < 2. 



L(t) 


n < o" 


n = o" 


n > a 


n + /i < 2 




(tint) 1 /* 


^l/(2+CT-n) 


n + = 2 


(t/lnt) 1 / 2 


t l/2 


(t/ln0 1/(2+<7 " n) 


n + /i > 2 


t l/(n+ M ) 




t l/(<7+^) 



Table 2: The growth law of the length scale L(t) for long-range forces, < a < 2, 
with various number of components, n, and conservation laws, /i. Note that n = d 
is excluded for n < a. 
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